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Making use of a derivative dispersion approach, we investigate the behavior of the the total 
cross section and the p parameter for pp and pp scattering from accelerator to cosmic ray energies. 
The discrepancies in the cosmic ray information is treated through the definition of two different 
ensembles of data. Simultaneous and individual fits to the above quantities through the Donnachie- 
Landshof and Kang-Nicolescu parametrizations allows to infer an upper bound for the intercept of 
the soft Pomeron and also show that the data investigated do not favor the Odderon hypothesis. 



I. INTRODUCTION 



o 
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^ , Forward quantities (total cross section, p, slope and curvature parameters) are fundamental quantities in the 
investigation of hadron-hadron scattering. In particular, total cross section and p are expressed in terms of the real 
and imaginary parts of the forward scattering amplitude by 
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' where, ^/s is the center of mass energy and t the four-momentum transfer squared. 

The theoretical description of these quantities as function of the energy constitute, presently, a great challenge. 
Among the variety of approaches, the analytical models are characterized by suitable analytic parametrizations to 
both quantities and the use of dispersion relations connecting real and imaginary parts of the amplitude (Principle of 
' Analyticity) . Also, particle-particle and particle-antiparticle scattering are simultaneously investigated by combining 
, even and odd amplitudes (Principle of Crossing Symmetry). 
O ■ Presently, the highest energies reached in accelerator experiments correspond to antiproton-proton scattering, ~ 
I ' 1.8 TeV and this value is nearly thirty times the highest energy reached in the case of proton-proton scattering, 
y/s = 62.5 GeV. These differences are the responsible for the difficulties present in the simultaneous investigation of 
' both reactions, meanly in the case of the behaviors as function of the energy and asymptotic limits. 
• • . In this communication we treat pp and pp forward scattering through analytical models and the central point is a 
^ ' critical investigation based on two observations: 

X 



1. The total cross section is an experimental quantity, basically determined by counting the number of elastic and 
' inelastic events in the scattering process. On the other hand, p is a free fit parameter, obtained through fits to 

the differential cross section data in the region of interference Coulomb - nuclear Therefore, the values 

of p extracted from experiments are model dependent. We conclude that atot and p have different status as 
physical quantities and this puts limitations in the simultaneous investigation of both quantities in constrained 
approaches, such as dispersion relations. 

2. Beyond the energy region of the accelerators there are experimental information on pp total cross section from 
cosmic ray experiments. Despite of the high error bars in the numerical results and also some discrepancies, 
they are usually shown as a support to several model and analysis predictions. 

The first observation suggests the investigation of the effect of simultaneous and independent fits to atot{s) and p{s) 
and the second, the determination of possible statistical bounds on the raise of these quantities, due to all the cosmic 
ray information presently available (taking account of the large error bars). These are the points we are interested 
in. The basic procedure is the choice of suitable parametrizations for atot with different asymptotic behaviors and 
the use of analyticity (dispersion) relations in order to connect cJtotis) and p{s). Observation 1 is treated through 
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individual and simultaneous fits to both quantities and observation 2 by selecting different ensembles of data based 
on distinct cosmic ray information, as will be explained in detail. The paper is organized as follows. In Sec. II we 
present the experimental data investigated and discuss the discrepancies concerning cosmic ray information. In Sec. 
Ill we review the main formulas in the dispersion relation approach and introduce two different parametrizations for 
the total cross section, namely the Donnachie-Landshof and the Kang-Nicolcscu models. The results of all the fits are 
presented in Sec. IV and the conclusions and some final remarks are the contents of Sec. V. 



II. EXPERIMENTAL DATA AND ENSEMBLES 



A. Accelerator Data 



In the case of accelerator experiments, data on atot from pp and pp total cross sections and also estimations of 
the p parameter have been accumulated for a long time. The data sets, compiled and analyzed by the Particle Data 
Group, have become a standard reference and are available (also readable files) at http://pdg.lhl.goi. Recent analysis 
has shown that general fits are stable for -^s above ~ 9 GeV . For this reason in this work we make use of these 
sets for energies above 10 GeV and in our analysis the statistic and systematic errors were added linearly. The total 
cross section data are shown in Fig. 1. 



B. Cosmic ray information 



The extraction of the proton-proton total cross section from cosmic ray experiments is based on the determination of 
the proton-air production cross section from analysis of extensive air showers. Detailed review on the subtle involved 
may be found in References [|j and Here we only recall the two main steps, stressing the model dependence 
involved. 

The first step concerns the determination of the proton-air production cross section, namely, the "inelastic cross 
section in which at least one new hadron is produced in addition to nuclear fragments" This is obtained by the 
formulas 

prod I , \ 2.4 X 10^ / I 2\ ^att 

("^^^ ^ A,_...(g/cm^) ' ^r>-aU9lcm ) = — , (2) 

where, Xp-air is the interaction length of protons in the atmosphere, Xatt the shower attenuation length and A: is a 
measure of the dissipation of the energy through the shower. The Xatt is a experimental quantity determined through 
the Xmax Attenuation Method or the Zenith Angle Attenuation Technique g]. On the other hand, the parameter k 
is obtained through Monte Carlo simulation and so is strongly model dependent 

In a second step cr*°* is obtained from o-p^/'^ir through the Multiple Diffraction Formalism by Glauber and Matthiae 
and taking account of the screening correction (Act): 

^p—air ^p—air ' ^ p—air ' p—air ' p—air ' 

Two important inputs at this point are the nucleon distribution function and the slope parameter _B, necessary in 
the parametrization of the proton-proton scattering amplitude: 

Fpp{s, q) cx exp{— ^^^}, (3) 

where is the four-momentum transfer squared. 

The results presently available from cosmic ray experiments cover the region ^/s ^6 — 40 TeV and are characterized 
by discrepancies, mainly due to the model dependence of k and B{s). However, we can classify these results according 
to the degree of consistence of the analysis and, simultaneously, by value of the total cross section extracted, as 
explained in what follows. 

From one hand the results usually quoted in the literature concern those obtained by the Fly's Eye Collaboration 
in 1984 and Akeno (Agasa) Collaboration in 1993 ||]. In the first case, the authors used k — 1.6 and a relation of 
proportionality between B and cr*°* as predicted by the Geometrical Scaling model ||l^. At ^/s — 30 TeV the cross 
sections results are 131: 



2 



prod 



540 ± 50 mb 



pp 



122 ± 11 mb 



The Akeno Collaboration used fc = 1.5 and the Durand and Pi method 12| to extract the proton-proton cross 
sections in the region ^/s ~ 6 — 20 TeV. The results are displayed in Fig. 1 together with the Fly's Eye value at 
= 30 TeV. 

On the other hand, in 1987 Gaisser, Sukhatme, Yodh (GSY) making use of the Fly's Eye result for (Jp^_^^^,. and the 
Chou - Yang prescription for B{s) obtained at ^/s = 40 TeV [Oj: 



■27 
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In 1993, Nikolaev claimed that the Akeno results should be corrected in order to take account of the differences 
between absorption and inelastic cross sections, leading to an increase of the published results by ~ 30 mb ||l^ . 

All these cosmic ray estimates together with the accelerator results are displayed in Fig 1 and show that, despite 
of the large error bars, we can identify two distinct set of estimations: one represented by the results of the Fly's Eye 
Collaboration together with those by the Akeno Collaboration; the other by the results of Gaisser, Sukhatme, Yodh 
with those by Nikolaev. Taken separately these two sets suggest different behaviors for the increase of the total cross 
section. 

At this point it should be stressed that the Fly's Eye Collaboration used the Geometrical Scaling hypothesis, which 
is violated even at the collider energy and so this result is certainly wrong. Durand and Pi asser ted in Ref. |l|] 
that their results presented in |l^], and used by the Akeno Collaboration, should be disregarded due to a wrong 
approximation. On the other hand, to our knowledge, the results by Gaisser, Sukhatme, Yodh and Nikolaev have 
never been critized in the literature. But we must yet take account of the model dependence involved as discussed in 

il- 

For the above reasons we shall quantitatively investigate the statistical determination of possible lower or upper 
limits in the increase of the total cross section, associated with the results by Akeno - Fly's Eye and Nikolaev - GSY, 
respectively. To this end we consider two ensembles of data, above 10 GeV, defined by: 

Ensemble I: pp accelerator data and pp accelerator data + Akeno + Fly's Eye. 
Ensemble II: pp accelerator data and pp accelerator data + Nikolaev + GSY. 
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FIG. 1. Total cross sections {pp and pp above -^s = 10 GeV): accelerator data and cosmic ray information available. 



III. ANALYTICAL APPROACH 



A. Analyticity relations 

For pp and pp scattering, crossing symmetry is expressed by 
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cr±(s) 
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and Analyticity allows to connect atot{s) and p(s) through two compact and symmetric formulas, 

raZ{s)^E{a+)~0{a^), 



(6) 
(7) 



where, E{fj^) and 0(ct_) are analytic transforms relating real and imaginary parts of Even (+) and Odd ) am- 
plitudes, respectively. They are generally represented by integral formulas (Integral Dispersion Relations) and in the 
case of one subtraction they are expressed by [p| 
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where K is the subtraction constant. 

At high energies we can take the limit sq — > and in this case they are equivalent to derivative forms given formally 

by [pi-ra 
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Operationally these expressions may be evaluated through the expansions |16| , |17| , |19[ : 
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With this approach, p(s) may be determined through suitable input parametrizations for the total cross sections. 
In this communication we shall use two different choices for these inputs: 
• Donnachie-Landshoff (DL) 

The Donnachie-Landshoff parametrization for the total cross sections are expressed by l20| 



(10) 



where, originally, X = 21.7 mb, Y = 56.08 mb, Z = 98.39 mb, e = 0.0808, r/ = 0.4525. 
For this parametrization Eqs. (8) and (9) reads 
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We see that this choice imphes 



Aa - aZ{s) - ctZ{s) = [Z - Y)s-^' ^ (asymptoticaUy) 

• Kang-Nicolescu (KN) 

The parametrizations for the total cross sections under the possibihty of a maximal Odderon contribution are given 

by [0 

aZis) = Ai+ Bilns + k In^ s, aZ{s) = A2 + In s + /c In^ s + 2Rs~^/'^ (11) 
and Eqs. (8) and (9) give 



aZ \ s 2 \ 2 J \ n J V 27r 



\ s 2 \ 2 J \ TT J V 27r 



Hot 

In this case we have 

Act = aZ{s) ~ aZtis) = {A2 - Ai) + {B2 - Bi) Ins + 2Rs-^f^ AA + ABlns 
so that if ^ and/or AS ^ then, asymptotically, Act 7^ 0. 

IV. FITS AND RESULTS 

As referred to before, we shall consider several cases: for both parametrizations (DL and KN) we perform indepen- 
dent and simultaneous fits to (Jtot{s) and p{s) by considering the two ensembles I and II, separately. 

In the individual procedure we first fit the atot data through Eqs. (10) (DL) or (11) (KN) and then fit the p data 
by letting free the subtraction constant K only; in the simultaneous procedure the subtraction constant X is a free 
parameter in the overall fit to both quantities. The fits were performed through the CERN-minuit routine. 

The results with the DL parametrization are shown in Figs. 2 and 3 and the statistical information are displayed 
in Table 1. Those with the KN parametrization in Figs. 4 and 5 and Table 2. 

V. CONCLUSIONS AND FINAL REMARKS 

From Tables 1 and 2 and Figures 2 to 5 we are lead to the following main conclusions: 

• Ensembles I and II 

In general, the results with ensembles I and II do not differ substantially and the cosmic ray information in ensemble 
II (Nikolaev and GSY) are not described in all the cases. Certainly this is due to the small number of points and 
the large error bars as compared with the accelerator data and also to the choices for the parametrizations. However 
noticeable differences occur with the Kang Nicolescu parametrization (Figures 4 and 5). 

• Simultaneous and individual fits 

The best statistical results are obtained through the simultaneous fits. For example, in the case of DL parametriza- 
tion, /d-o.f. ~ 1.0 (simultaneous) and ^ 1.7 (individual for p). This is expected since in the simultaneous case we 
have more statistical information (degree of freedom). 

On the other hand simultaneous fits constraint the rise of the total cross section. For example, in the case of the 
DL parametrization with ensembles I and II we obtained: 

ef„rf, ~ 0.088 ^ ei^^^i ~ 0.083 (reduction ^ 6%) 

- 0-091 - efLn; ~ 0.084 (reduction ^ 8%) 

As commented before, atot and p do not have the same status as physical quantities since the p value is model 
dependent. Therefore we understand that the simultaneous fit may led to erroneous predictions for atot in the 
asymptotic region (underestimation). 
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• Upper bounds for the soft Pomeron intercept 

The upper bound was obtained with ensemble II in the case of individual fits 

l + e = 1.091± 0.003 1 - 77 = 0.65 ±0.03 

• Odderon 

From Tables 3 and 4, and AB ^ in 3 cases and the only exception is for individual fits with ensemble II, for 
which AA = 5 ± 3 and AB ^ 0. This predicts a crossing with crf^j becoming higher than af^^. However, in this case 
p is not described as shown in Fig. 5. We conclude that our results do not favor the Odderon possibility. 

We observe that a similar and interesting analysis, taking quantitative account of the cosmic ray information, has 
been presented in this Workshop by Carlos Avila . Although the technical approach is different the global results 
are in agreement with those presented here and in Ref . ||l9| . 
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TABLE I. Simultaneous and individual fits to atot and p for Ensembles I 

and II with the DL paraiiictrizatioii. 



Fit: 


Simultaneous 


Individual 


Ensemble: 


I 


II 


I 


II 


N. points 


165 


165 


102 


102 


x'/d-o.f. 


0.84 


0.98 


0.76 


0.96 


X (ml)) 


21.4 ±0.4 


21.1 ±0.4 


20.0 ± 0.7 


19.3 ± 0.7 


1 (mbj 


()7 = () 


()7 ^ :> 


18 ^ 


I() ^ :> 


Z (mb) 


114± 11 


112 ± 11 


75± 10 


70 ±9 


V 


0.48 ± 0.02 


0.47 ±0.02 


0.37 ±0.04 


0.35 ± 0.03 


e 


0.083 ± 0.002 


0.084 ± 0.002 


0.088 ± 0.003 


0.091 ± 0.003 


K 


306 ± 54 


307 ± 54 






N. points (p) 






63 


63 


x'/d.o.f. 






1.55 


1.84 


K 






235 ± 32 


245 ± 32 



TABLE II. Simultaneous and individual fits to atot and p for Ensembles I 
and II with the KN parametrization. 



Fit: 


Simultaneous 


Individual 


Ensemble: 


I 


II 


I 


II 


N. points 


165 


165 


102 


102 


x'/d.o.f. 


0.77 


0.86 


0.78 


0.91 


Ai (mb) 


44.4 ± 0.7 


45.0 ±0.7 


45 ± 1 


47 ± 1 


A2 (mb) 


44.5 ±0.7 


45.1 ±0.7 


45 ±3 


52 ±3 


Bi (mb) 


-2.9 ±0.2 


-3.1 ±0.2 


-2.9 ±0.3 


-3.6 ±0.4 


B2 (mb) 


-2.9 ±0.2 


-3.1 ±0.2 


-2.9 ±0.6 


-4.2 ±0.6 


k (mb) 


0.33 ±0.01 


0.35 ±0.01 


0.33 ±0.06 


0.39 ±0.03 


R (mb GeV) 


25.3 ±0.9 


25.4 ±0.9 


24 ±7 


12 ±7 


K 


-72 ± 46 


-63 ± 46 






N. points (p) 






63 


63 


x'/d.o.f. 






1.05 


35.2 


K 






-198 ±32 


-1369 ± 32 
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FIG. 2. Fits to atot with the DL parametrization: pp (solid) andpp (dashed) and ensembles I (blue) and II (red). Simultaneous 
(up) and individual (down). 
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FIG. 3. Fit to p with the DL parametrization: pp (solid) and pp (dashed) and ensembles I (blue) and II (red). Simultaneous 
(up) and individual (down). 
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FIG. 4. Fits to atot with the KN parametrization: pp (solid) and pp (dashed) and ensembles I (blue) and II (red). Simulta- 
neous (up) and individual (down). 
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FIG. 5. Fit to p with the KN parametrization: pp (solid) and pp (dashed) and ensembles I (blue) and II (red). Simultaneous 
(up) and individual (down). 
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